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Abstract 

We analyze the issue of "particle motion" in general relativity in a systematic and rigorous 
way by considering a one-parameter family of metrics corresponding to having a body (or black 
hole) that is "scaled down" to zero size and mass in an appropriate manner. We prove that the 
limiting worldline of such a one-parameter family must be a geodesic of the background metric and 
obtain the leading order perturbative corrections, which include gravitational self-force, spin force, 
and geodesic deviation effects. The status the MiSaTaQuWa equation is explained as a candidate 
"self-consistent perturbative equation" associated with our rigorous perturbative result. 
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It is of considerable interest to determine the motion of a body in general relativity in 
the limit of small size, taking into account the deviations from geodesic motion arising from 
gravitational self-force effects. There is a general consensus that the gravitational self-force 
is given by the "MiSaTaQuWa equations" : In the absence of incoming radiation, the motion 
is given by 

«"V*u" = -\isT + uV)(2VXf ~ V»h%% {T) u<>u° , (1) 

= M (<WV - ^ + /mv) (x, z{t')) u»'u»'dT> , (2) 

where G + is the retarded Green's function for the wave operator V a V a /i Ml , — 2R a flu /3 h a /3. 
(Note that the r~ limit of integration indicates that only the part of G + interior to the light 
cone contributes to h^.) However, all derivations contain some unsatisfactory features. 
This is not surprising in view of the fact that, as noted in "point particles" do not make 
sense in nonlinear theories like general relativity! 

• Derivations that treat the body as a point particle require unjustified "regularizations" . 

• Derivations using matched asymptotic expansions Q make a number of ad hoc and/or 
unjustified assumptions. 

r-| 

• The axioms of the Quinn-Wald axiomatic approach [3| have not been shown to follow 
from Einstein's equation. 

• All of the above derivations employ at some stage a "phoney" version of the linearized 
Einstein equation with a point particle source, wherein the Lorenz gauge version of 
the linearized Einstein equation is written down, but the Lorenz gauge condition is 
not imposed. 

How should gravitational self-force be rigorously derived? A precise formula for gravita- 
tional self-force can hold only in a limit where the size, R, of the body goes to zero. Since 
"point-particles" do not make sense in general relativity — collapse to a black hole would 
occur before a point-particle limit could be taken — the mass, M, of the body must also go 
to zero as R — > 0. In the limit as R, M — > 0, the worldtube of the body should approach a 
curve, 7, which should be a geodesic of the "background metric". The self-force should arise 
as the lowest order in M correction to 7. In the following, we shall describe an approach 
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that we have recently taken to derive gravitational self-force in this manner. Details can be 
found in (4]. 

The discussion above suggests that we consider a one-parameter family of solutions to 
Einstein's equation, (g^(A), T /U/ (A)), with R(X) — > and M(A) — ► as A — > 0. But, what 
conditions should be imposed on (g^ v (X), T Ml/ (A)) to ensure that it corresponds to a body that 
is shrinking down to zero size, but is not undergoing wild oscillations, drastically changing 
its shape, or doing other crazy things as it does so? 

As a very simple, explicit example of the kind of one-parameter family we seek, consider 
the Schwarzschild-deSitter metrics with M = A, 

ds ^X) = -(1- — - Cr 2 )dt 2 + (1 - — - Cr 2 Y l dr 2 + rW. (3) 

If we take the limit as A — > at fixed coordinates (t, r, 9, <fi) with r > 0, it is easily seen that 
we obtain the deSitter metric — with the deSitter spacetime worldline 7 defined by r = 
corresponding to the location of the black hole "before it disappeared" . However, there is 
also another limit that can be taken. At each time to, one can "blow up" the metric <7 M!/ (A) 
by multiplying it by A -2 , i.e., define 

V(A) = \-\ u (\). (4) 

We correspondingly rescale the coordinates by defining f = r / A, t = (t — t )/X. Then 

ds 2 (X) = -(1 - 2/f - X 2 Cf 2 )dt 2 + (1 - 2/f - X 2 Cf 2 )- 1 df 2 + fW (5) 

In the limit as A — > (at fixed (t,r,9,<f>)) the "deSitter background" becomes irrelevant. 
The limiting metric is simply the Schwarzschild metric of unit mass. The fact that the limit 
as A — >• exists can be attributed to the fact that the Schwarzschild black hole is shrinking 
to zero in a manner where, in essence, nothing changes except the overall scale. 

The simultaneous existence of both of the above types of limits charaterizes the type of 
one-parameter family of spacetimes g^iX) that we wish to consider. More precisely, we wish 
to consider a one parameter family of solutions 57^ (A) satisfying the following properties: 

• (i) Existence of the "ordinary limit": There exist coordinates x a such that g^X, x a ) is 
jointly smooth in (A, x a ), at least for r > RX for some constant R, where r = a/XX x *) 2 
(i = 1,2,3). For all A and for r > .RA, g^ v {X) is a vacuum solution of Einstein's 
equation. Furthermore, g^ v {X = 0,x a ) is smooth in x a , including at r = 0, and, for 
A = 0, the curve 7 defined by r = is timelike. 
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• (ii) Existence of the "scaled limit": For each to, we define t = (t — to) A, x % = x l /\. 
Then the metric g^p(A; to; x a ) = A~ 2 g^p(A; t ; x") is jointly smooth in (A,t ;x Q ) for 
f = r/A > R. 

The above two conditions must be supplemented by an additional "uniformity require- 
ment", which can be explained as follows. From the definitions of gp& and x^, we can relate 
coordinate components of the barred metric in barred coordinates to coordinate components 
of the unbarred metric in corresponding unbarred coordinates, 

<7pp(A; t ; t, x l ) = g^(X; t + At, Xx l ) . (6) 

Now introduce new variables a = r and (5 = A/r = 1/f, and view the metric components 
g^u(X) as functions of (a, (5, t, 9, 0), where 9 and are defined in terms of x % by the usual 
formula for spherical polar angles. We have 

gp,„(a(3,t ; i, 1/(5, 9,(j>) = g^{a(5,t = t + At; a,9,(b) . (7) 

Then, by assumption (ii) we see that for < (5 < 1/R, g^ u is smooth in (a, (3) for all a 
including a = 0. By assumption (i), we see that for all a > 0, g^ u is smooth in (a, (5) for 
(5 < 1/R, including (5 = 0. Furthermore, for (5 = 0, g^ v is smooth in a, including a = 0. 

We now impose the additional uniformity requirement on our one-parameter family of 
spacetimes: 

• (iii) g^y is jointly smooth in (a, (5) at (0,0). 

We already know from our previous assumptions that g^ (A; to, r, 9, (ft) and its derivatives 
with respect to x a approach a limit if we let A —>■ at fixed r and then let r — > 0. The 
uniformity requirement implies that the same limits are attained whenever A and r both go 
to zero in any way such that A/r goes to zero. 

It has recently been proven in 5j that an analog of the uniformity requirement holds for 
electromagnetism in Minkowski spacetime in the following sense: Consider a one-parameter 
family of charge-current sources of the form J 71 (A, t, x % ) = J^(X, t, x l / A) where is a smooth 
function of its arguments and x % = defines a timelike worldline. Then the retarded solution, 
F^iX, x^), is a smooth function the variables (a, (5, t, 9, (b) in a neigborhood of (a, (5) = (0, 0). 
In the gravitational case, we do not have a simple relationship between the metric and the 
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stress-energy source, and in the nonlinear regime, it would not make sense to formulate the 
uniformity condition in terms of the behavior of the stress-energy. Consequently, we have 
formulated this condition in terms of the behavior of the metric itself. 

The uniformity requirement implies that the metric components can be approximated 
near (a, (3) = (0, 0) with a finite Taylor series in a and j3, 

^(A;t,r,0,0) = ^)^r n f-J {a^) nm {t,e, 0), (8) 

n=0 m=0 

where remainder terms have been dropped. This gives a far zone expansion. Equivalently, 
we have 

N M / l\ m 

9pp(X;t ;i,f,6,(f>) = ^^(Af) n f-J (<v) nro (t + Xt,6, <j>) . (9) 

n=0 m=0 ^ ' 

Further Taylor expanding this formula with respect to the time variable yields a near zone 
expansion. Note that since we can express g^ at A = as a series in 1 jf as f — > oo and since 
g^ v at A = does not depend on i, we see that <^p(A = 0) is a stationary, asymptotically 
flat spacetime. 

The curve 7 to which our body shrinks as A —>■ (see condition (i) above) can now be 
proven to be a geodesic of the metric g flv (X = 0) as follows: Choose the coordinates x a 
so that at A = they correspond to Fermi normal coordinates about the worldline 7. In 
particular, we have g^ v = rj^ on 7 at A = 0. It follows from ([8]) that near 7 (i.e., for small 
r) the metric g^ v must take the form 

g»„ = V,u + 0{r) + A f ^hM + 0(1)) + 0(A 2 ) (10) 

Now, for r > 0, the coefficient of A, namely 

V = ^+o(i) (ii) 

must satisfy the vacuum linearized Einstein equation off of the background spacetime 
9ixu(X = 0). However, since each component of h^ v is a locally L l function, it follows 
immediately that h^ u is well defined as a distribution. It is not difficult to show that, 
as a distribution, satisfies the linearized Einstein equation with source of the form 
N^S^ix*), where is given by a formula involving the limit as r — ► of the angular 
average of C^ u and its first derivative. The linearized Bianchi identity then immediately 
implies that N^ u is of the form Mu^u u with M constant, and that 7 is a geodesic for 



Our main interest, however, is not to rederive geodesic motion but to find the leading 
order corrections to geodesic motion that arise from finite mass and finite size effects. To 
define these corrections, we need to have a notion of the "location" of the body to first order 
in A. This can be defined as follows: Since g^p(A = 0) is an asymptotically flat spacetime, its 
mass dipole moment can be set to zero (at all to) as a g au g e condition on the coordinates x l . 
The new coordinates x l then have the interpretation of being "center of mass coordinates" 
for the spacetime g^p(A = 0). In terms of our original coordinates x a , the transformation to 
center of mass coordinates at all t corresponds to a coordinate transformation x a — > x a of 
the form 

x a (t) = x a -XA a {t,x i ) + 0(X 2 ). (12) 

To first order in A, the world line defined by x l = should correspond to the "position" 
of the body The first-order displacement from 7 in the original coordinates is then given 
simply by 

Z\t) = A l (t,x j =0). (13) 

The quantity Z % is most naturally interpreted as a "deviation vector field" defined on 7. 
Our goal is to derive relations (if any) that hold for Z % that are independent of the choice 
of one-parameter family satisfying our assumptions. 

We now choose the x a coordinates — previously chosen to agree with Fermi normal coor- 
dinates on 7 at A = — to correspond to the Lorenz/harmonic gauge to first order in A. To 
order A 2 , the leading order in r terms in g a p are, 

g a p(\; t, x*) = 7] aP + B a0j (t)x i x j + 0(r 3 ) 

+ A (j^Sap + h$(t, 0) + h^(t, OK + MTl aP {t) + 0(r 2 )^j 

( M 2 2 1 

+ A 2 ( — (-2t a t p + 3n a np) + -P t (t)n l 5 af3 + -t (a S p)j {t)n' 

+ l -K a p{t, 9, 0) + H aP {t, 9, 0) + 0(r)^j + 0(A 3 ) 

Here B and TZ are expressions involving the curvature of 57^ (A = 0) and we have introduced 
the "unknown" tensors K and H. The quantities P l and S a p turn out to be the mass 
dipole and spin of the "near- zone" background spacetime ^p(A = 0). For simplicity, we 
have assumed no "incoming radiation". Hadamard expansion techniques and 2nd order 
perturbation theory were used to derive this expression. 



(14) 
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Using the coordinate shift x A 



to cancel the mass dipole term, the above 



expression translates into the following expression for the scaled metric 

2M M 2 1/1 
9ap(to) = Vap + —8 a p + — {-2t a t p + 3n a np) + —t^S^n 1 +01 — 



+ A 
+ A 2 



Kl + 2A^) + + I? VV*'' + o(±)+io ( I 



(a,/3)7^' 



+ + -K aP + —t^S^n 3 + O 



0(A 3 



(15) 



The terms that are first order in A in this equation satisfy the linearized vacuum Einstein 
equation about the background "near zone" metric (i.e., the terms that are Oth order in 
A). From this equation, we find that dSij/dt = 0, i.e., to lowest order, spin is parallelly 
propagated along 7. 

The terms that are second order in A in this equation satisfy the linearized Einstein 
equation about the background "near zone" metric with source given by the second order 
Einstein tensor of the first order terms. Extracting the £ = 1, electric parity, even- under- 
time-reversal part of this equation at 0(l/f 2 ) and 0(i/f 3 ), we obtain (after considerable 
algebra!) 

z i,oo = 7^jS kl R k mi - RojOiZ j - (h}^ - -/loo)* J • (16) 
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In other words, in the Lorenz gauge, the deviation vector field, Z a , on 7 that describes the 
first order perturbation to the motion satisfies 



u c V c {u b V b Z a ) 



2M' 



Rbcd a S bC U d ■ 



R bc /u b u d Z c -(g ab + u a u b )(V d htf-^V b hlf)u c u d . (17) 



Equation ffTTl) gives the desired leading order corrections to motion along the geodesic 7. 
The first term on the right side of this equation is the Papapetrou "spin force" , which is the 
leading order "finite size" correction. The second term is just the usual right hand side of the 
geodesic deviation equation; it is not a correction to geodesic motion but rather allows for 
the possibility that the perturbation may make the body move along a different geodesic. 
Finally, the last term describes the gravitational self-force that we had sought to obtain, 
i.e., the corrections to the motion caused by the body's self- field. Equation (fT7j) gives the 
correct description of motion when the metric perturbation is in the Lorenz gauge. When 
the metric perturbation is expressed in a different gauge, the force will be different {4]. 



Although we have now obtained the perturbative correction to geodesic motion due to 
spin and self-force effects, at late times the small corrections due to self-force effects should 
accumulate (e.g., during an inspiral), and eventually the orbit should deviate significantly 
from the original, unperturbed geodesic 7. When this happens, it is clear our perturbative 
description in terms of a deviation vector defined on 7 will not be accurate. Clearly, going to 
any (finite) higher order in perturbation theory will not help (much). However, if the mass 
and size of the body are sufficiently small, we expect that its motion is well described locally 
as a small perturbation of some geodesic. Therefore, one should obtain a good description 
of the motion by making up (!) a "self- consistent perturbative equation" that satisfies the 
following criteria: (1) It has a well posed initial value formulation. (2) It has the same 
"number of degrees of freedom" as the original system. (3) Its solutions correspond closely 
to the solutions of the of the original perturbation equation over a time interval where the 
perturbation remains small. In some sense, such a self-consistent perturbative equation 
would take into account the important ("secular") higher order perturbative effects (to all 
orders), but ignore other higher order corrections. Such equations are commonly considered 
in physics. The MiSaTaQuWa equations appear to be a good candidate for a self-consistent 
perturbative equation associated with our perturbative result. 

In summary, we have analyzed the motion of a small body or black hole in general 
relativity, assuming only the existence of a one-parameter family of solutions satisfying 
assumptions (i), (ii), and (iii) above. We showed that at lowest ("zeroth") order, the motion 
of a "small" body is described by a geodesic, 7, of the "background" spaceetime. We then 
derived is a formula for the first order deviation of the "center of mass" worldline of the body 
from 7. The MiSaTaQuWa equations then arise as (candidate) "self-consistent perturbative 
equations" based on our first order perturbative result. Note that it is only at this stage 
that "phoney" linearized Einstein equations come into play. 

We have recently applied this basic approach to the derivation of self-force in electro- 
magnetism [jj , and have argued that the reduced order form of the Abraham-Lorentz-Dirac 
equation provides an appropriate self-consistent perturbative equation associated with our 
first order perturbative result (whereas the original Abraham-Lorentz-Dirac equation is ex- 
cluded). It should be possible to use this formalism to take higher order corrections to the 
motion into account in a systematic way in both the gravitational and electromagnetic cases. 
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